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Abstract 

The decay ^-'(25') rjj is used to measure, for the first time, all prominent ry-meson branching 
fractions with the same experiment in the same dataset, thereby providing a consistent treatment 
of systematics across branching fractions. We present results for 77 decays to 77, 7r+7r~7r°, 37r°, 
7r+7r~7, and e^e~^, accounting for 99.9% of all 77 decays. The precisions for several of the branching 
fractions and their ratios are improved. Two channels, 'k~^tt~^ and e^e^7, show results that differ 
at the level of three standard deviations from those previously determined. 
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The r] meson was discovered almost half a century ago It is the second-lightest 
meson, considered to consist of u, d, and s quarks, and studying its decays into pions, 
electrons, and photons gives insight into different aspects of non-perturbative QCD and 
electromagnetic phenomena. Measurements of the r] decay properties come from many 
different experiments, and almost all exclusive branching fraction determinations are made 
relative to other t] decays. The Particle Data Group (PDG) [3] uses 43 such measurements 
in a fit to determine the branching fractions to 77, 37r'^, 7r+7r~7r°, 7r+7r~7, 7r°77, e~^e~'y, 
fi'^fi~'y, and 7i^7i~e^e~ , as well as the total width. 

The analysis presented here studies rj decays in the reaction e~^e~ —>■ 'ip{2S) riJ/ip 
with 1] '-f'-f, 37r°, 7r+7r~7r°, 7r+7r^7 and e~^e~'~f. We identify the J/ip through its decays to 
e~^e~ and The choice of modes addresses the known branching fractions of (9(0.1%) 

and larger, and covers 99.88% of the rj decay modes when using the branching fractions 
from Ref. [2|]. The strength of this analysis lies in the simultaneous and similar treatment 
of charged and neutral r] decay products, cross-feed of different modes into each other, and, 
with the same analysis procedure, estimates of backgrounds from other XJ/ip sources. 

The CLEO-c detector is described in detail elsewhere [3]. Its features exploited here are 
the 93% solid angle coverage of precision charged particle tracking and an electromagnetic 
calorimeter consisting of 7784 CsI(Tl) crystals, the barrel portion of which has a vertex- 
pointing geometry. The barrel calorimeter and two open-cell drift chambers are concentric 
with the colliding beams and embedded inside a 1 T axial magnetic field provided by a 
superconducting solenoid. The small inner chamber has six cylindrical stereo layers (drift 
cells canted at an angle to the chamber axis), and the outer, larger chamber has 47 layers, 
the inner 16 of which are axial and the outer 31 stereo. (About 5% of the data used here 
were acquired in the earlier CLEO III detector configuration, which differed from CLEO-c 
primarily by having a four-layer silicon strip vertex detector in place of the inner tracking 
chamber.) The tracking system enables momentum measurements for particles with mo- 
mentum transverse to the beam exceeding 50 MeV/c, and achieves resolution ap/p ~0.6% 
at p=l GeV/c. The barrel calorimeter reliably measures photon shower energies down to 
E^=30 MeV and has a resolution of aE/E ~5% at 100 MeV and 2.2% at 1 GeV. 

The data sample comprises about 27M ip{2S) decays, corresponding to about O.IM rj 
decays produced with a J/ip ^ tag. 

We determine the detection efficiency and background levels with Monte Carlo (MC) 
samples that were generated using the EvtGen event generator [ol] and a GEANT-based [3| 
detector simulation. We model — > 77 and Sn^ according to phase space. The mode n~^7r~^ 
is simulated as mediated by a p° ^ tt+tt" decay, weighted with a factor ~ E^, where E^ 
is the photon energy in the f] center-of-mass system. We generate tt^it'tt^ according to the 
distribution measured in The simulation of e^e~7 is analogous to tt^ e+e^7 ("Dalitz 
decay") 0. 

The event selection proceeds as follows. We select the J/ip ^ i^i" track candidates 
within polar angles |cos6'^±| < 0.83, adding bremsstrahlung photons within a cone of 
lOOmrad around the track momentum vector at the collision point. We identify leptons 
through the ratio of energy deposition in the calorimeter associated with the track, E, to 
the track momentum measured in the drift chamber, p: For electron (muon) candidates, we 
require E/p values of > 0.85 (< 0.25) for one lepton and > 0.50 (< 0.50) for the other. 
We impose kinematic constraints by fitting the two lepton candidates to a common orig- 
inating vertex (where the figure of merit is given by Xj^^ ^/d.o.f.) and to the J/ip mass 
(Xj/^^m/d.o.f.). We keep candidates that have Xj/xii^v/^-'^-^- < 20 and Xj/^ m/'^-O-f- < 20, 
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which keeps signal decays with high efficiency, as evident from Fig. [H The photons in the 
f] decay products are required to be in the region of best calorimeter performance and least 
material in front of the crystals, | cos^^| < 0.75, and not be matched or close to a track's 
projection into the calorimeter. 

We then proceed to use kinematic constraints once more for improved event cleanliness: 



The fitted J/ip and the rj decay products are constrained, together with the beam spot [10 



to a common vertex, and then to the ip{2S) mass. This results in a very clean separation of 
final states. We apply mode-dependent restrictions on the quality of these fits, denoted by 
X^{2S),v/^-'^-^- Xv)(2S),m/'^-0-f-' respectively. Conversion events originating from — > 77 
decay can fulfil the e"*"e~7 pre-selection, but have a poor ip{2S) vertex fit; hence we apply 
a stricter cut of X^{2S),v/^-'^-^- < 4 in this channel (see Fig. [1]). All other modes require 
X^(25) v/'i-O-f- < 20. The mass fit has mode-dependent cuts, set as loosely as possible while 
preserving sample cleanliness: X^{25),m/*i-0-f- < 20 for tt+tt^tt^, < 10 for 77, Stt", and e"'"e~7, 
< 5 for 7r"'"7r~7. 

After this step, we define the following signal windows: p{J/ip) = 170 — 230MeV/c, and 
two ranges for m{r]): 542 — 554 MeV for 7r+7r~7, 535 — 560 MeV for all others. Final state 
specific characteristics are: (1) 77: > 200 MeV, to suppress photons from 4'{2S) 7Xci- 
(2) 7r+7r~7r°: We search for two photons with > 30 MeV and ^(77) = 100 — 160 MeV, and 
constrain them to the tc^ mass. (3) Stt^: We search for six photons, but do not attempt to 
make assignments to 7r° candidates because doing so typically results in multiple comparably 
probable assignments. (4) tt^tc"'-/: E^ > 100 MeV and m(7r"*"7r~) > 300 MeV. (5) e~^e~'-f: 
We add bremsstrahlung photons to the soft electrons as with J/ip —>■ and the soft 

tracks must satisfy |cos6'| < 0.8. In addition, we require m(e"^e~) < 300 MeV. A substan- 
tial number of 77 events with a conversion survive the vertex restriction described above 
and fake the e"'"e"7 signature; indeed, this type of background has necessitated substantial 
subtractions in previous measurements of this mode 0, 0]. These conversions tend to occur 
at the discrete locations such as the beam pipe and tracking chamber boundaries, but are 
reconstructed as if they originated at the interaction point. Consequently they create an 
artificial mass peak near 10 MeV as seen in the lower right of Fig. [2J We remove the mass 
region m{e^e~) = 8 — 20 MeV to suppress this background and the systematic uncertainties 
associated with it. 

For all five rj decay channels, we keep the two J /ip decay modes separate. The fit quality 
for data and simulation is compared in Fig. [H As a cross-check, we also perform the analysis 
without the ip{2S) kinematic fit: consistent results are obtained, but in most modes with 
far worse background conditions and larger uncertainties. 

The main backgrounds arise from cross-feed between the rj modes and from other ip{2S) — >■ 
XJ/ip transitions, mostly X = vr+vr", 7r°7r°, and 77 through Xcj- We select such exclusive 
event samples using selections similar to the rj signal decays, including the kinematic fits. 
Backgrounds from these XJ/tp channels into the t] signals are then determined by scaling 
the MC predictions so as to match the observed XJ/ip yields in data, and subtracted. 
The statistical uncertainties of these subtractions are accounted for. We find that 77, 37r°, 
7r+7r^7r°, and 7r+7r^7 have such backgrounds at the levels of 1-2%. Examination of the r] 
mass sidebands revealed no discrepancy between data yields and MC estimate; the only 
exception is 7r+7r~7, where data exceeds MC by an amount which, when extrapolated into 
the signal region, corresponds to a background of (2.8 ± 1.1)% and is subtracted in addition 
to the other Monte Carlo predictions. The mode e"*"e~7 has a background of about 5% 
due mostly to 77 conversions which survived the m{e~^e~) and tight vertex fit restrictions. 
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Other (non-J/'?/') il^{2S) decays do not fake the signal signature at any appreciable level. We 
use a 20.7 pb~^ sample of data taken at a center-of-mass energy of 3.670 GeV to estimate 
continuum background (scaled by luminosity and energy dependence), which is found to be 
negligible. 

All inspected experimental observables show good agreement between data and the sum 
of our MC samples, normalized according to their relative population in the data. A selection 
of comparisons is shown in Figs. [2] and [31 

Table [I lists observed yields and the estimated background. We observe significant, clean, 
and unambiguous signals for all our target modes. 

Our measurements are performed as ratios between efficiency-corrected event yields of 
pairs of r] final states, separately for J/tp —>■ e~^e~ and ;U^/i~. This allows cancellation of 
all lepton-related systematic uncertainties, such as track finding, lepton identification, and 
J/ip fitting. We then proceed to combine the two measurements of each ratio, where the 
ry-related uncertainties are treated as fully correlated. We note that the absolute detection 
efficiency for iIj{2S) — > rjj/ilj, r] — > 77, J/ip ^ is about one third. 

Sources of systematic uncertainty and the values assigned are: Track finding (0.3% per 
track, added linearly photon finding (0.4% per photon, added linearly |11]), sideband 
subtraction (1.1%, 7i~^7i~'y only), trigger (0.1-0.5%, mode-dependent), and MC statistics 
(0.4-1.0%, mode- dependent), other effects in the detector simulation (0.5%). We also make 
reasonable variations in decay modeling at the MC generator level and assign uncertainties 
accordingly: 0.1% for 37r° to account for the slight deviation from phase-space-prescribed 
decay observed in Ref. [l2|, 0.9% for tt+tt^tt^ based on the experimental uncertainty of 
the slope parameter in Ref. jsl, 3% for 7r+7r^7 to include a slightly different lineshape 
parameterization of the intermediate p° meson, and 5% for e^e~7 to allow for changes in 
the polar angle distribution of the from the rj decay and in the m{e~^e~) spectrum that 
remain consistent with our measurements. All uncertainties are added in quadrature, except 
where correlations between modes have to be observed. 

The results for ratios of branching fractions are shown in Table Ull The fo^' the ratios 
relative to 77 to agree between J/ip —>■ e~^e~ and fi~^fi~ is 5.9 for four degrees of freedom, 
corresponding to a confidence level of ~ 20%. We designate the following four ratios as 
constituting a complete set, having minimal systematic correlation with each other: Sn^/'j'j, 
7r+7r~7r°/77, 7r+7r^7/7r+7r^7r°, and e+e~7/7r+7r~7. We compare to the single most precise 
other measurements in Fig. |H 

Under the assumption that our five signal modes account for all of the t] decay modes, 
we combine the ratios between them to form absolute branching fraction measurements. 
Correlations between uncertainties are taken into account. Other possible rj decay modes 
are either forbidden and/or have been found to be below 0.2% in branching fraction 0]: 
We include 0.3% as a systematic uncertainty in the absolute branching fraction results. The 
results are presented in Table [Till together with those from PDG 2006 [2] for the global fit to 
all measurements. In all five modes, the statistical uncertainty is larger than or comparable 
to the systematic error. A visual comparison can be found in Fig. m 

To summarize, we have studied five rj decay modes using the decay chain i'{2S) — >• rjj/ip, 
J/ijj —>■ ee and fifi: rj — *• 77, 37r°, 7r+7r~7r°, 7r+7r~7, and e+e~7. We have presented ratios 
between these modes as well as absolute rj branching fractions to these final states. This 
is the first analysis that covers this range of 1] decay modes, summing up to 99.9% of the 
known rj decays, and determines their absolute branching fractions in the same experiment. 
Several of the relative and derived absolute branching fractions obtained in this analysis are 
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either the most precise to date or first measurements. In particular, we note that our result 
for 7r+7r~7 is about ~ 15% (3.2cr) smaller than previous measurements, and for e^e~7 is 
~ 57% (2.9cr) larger. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 
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TABLE L For each rj decay channel, the observed yields in the ip{2S) on-resonance sample (A^^^S)^^ 
background from cross-feed between r/ modes {N'^^), and background from other XJ/ijj decays 
(^jyXj/ip^ separately for J/ip ^ ee and J/V' fJ'l-i'- 



Channel 










J/ip e+e 




J 1^1, ^ ^l+ 




77 


6324 - - 


66 


7376 - - 


114 




850 - - 


18 


1004 - - 


15 


vr+TT^vr'' 


1884 - 4 - 


12 


2052 - 5 - 







403 - 3 - 


17 


498 - 2 - 


20 


e+e~7 


82 - 4 - 





100 - 6 - 






TABLE IL Ratios of r] branching fractions. For each combination, the efficiency ratio, separately 
for J — e+e~ and J /il) —i- fJ^^lJ-' , the level of consistency between the J/ip —>■ e+e~ and fj.'^n~ 
result, expressed in units of Gaussian standard deviations, cr^^/^^, and the combined result for the 
branching ratio. The dagger symbol indicates that this result is most precise measurement to date. 



Channel 


eff. ratio 




branching fraction ratio 






ee 






371*^/77 


0.15 


0.15 


1.0 


0.884 ± 0.022 ± 0.019 


7r+7r^7r*'/77 


0.50 


0.49 


-2.2 


0.587 ± 0.011 ± 0.009t 


7r+7r^7/77 


0.63 


0.60 


0.2 


0.103 ± 0.004 ± 0.004t 


e+e~7/77 


0.53 


0.52 


0.1 


0.024 ± 0.002 ± O.OOlt 


37r°/7r+7r-7r° 


0.30 


0.32 


2.1 


1.496 ± 0.043 ± 0.035t 


vr + vr ~ 7/7r+ vr ^ vr'^ 


1.27 


1.24 


1.1 


0.175 ± 0.007 ± 0.006 


e+e~7/7r+7r~7r'^ 


1.07 


1.06 


0.5 


0.041 ± 0.003 ± 0.002t 


e+e~7/7r+7r~7 


0.84 


0.86 


0.0 


0.237 ± 0.021 ± 0.015 



TABLE III: For each r] decay channel, absolute branching fraction measurements for J /ip e+e~ 
and J/^lJ l^'^l^'^ combined, with statistical and systematic uncertainties (middle column), as 
determined in this work. The last column shows the PDG fit result Q]- All but 77 are first 



measurements 
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Channel 


this work (%) 


PDG [2] (%) 


77 


38.45 ± 0.40 ± 0.36 


39.38 ± 0.26 


3^0 


34.03 ± 0.56 ± 0.49 


32.51 ± 0.28 


vr+vr^Tr" 


22.60 ± 0.35 ± 0.29 


22.7 ± 0.4 


7r+7r~7 


3.96 ± 0.14 ± 0.14 


4.69 ± 0.11 


e+e~7 


0.94 ± 0.07 ± 0.05 


0.60 ± 0.08 
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FIG. 1: Top to bottom: 77, 3tt^, tt'^'tt^tt'^, 7r"*"7r~7, e'^e~7. For each channel, left to right: 
Goodness-of-fit for J/tp vertex and mass fit, and for tp{2S) vertex and mass fit. Points: data. 
Dotted line: Signal MC. Solid line: Sum of all MC. Arrows indicate selection requirements. Cuts 
have been applied to all quantities with the exception of the one plotted. 
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2540607-002 




m(7r+7r-)(GeV) Ey(GeV) 




m(e+e-)(GeV) logio(m(e+e-) (MeV) ) 

FIG. 2: Distributions for individual channels. Top row: Svr^, highest and lowest photon energies. 
Second row: 'K~^7:~7r^, kinematic distribution of the three pions, and two-photon invariant mass; y 
is a function of the kinetic energy of the vr*^ (Tq) and the sum of the kinetic energies of all pions 
(Q)- y = (37o/Q) — 1. Third row: 7r"'"7r~7, invariant mass of the two pions, and photon energy. 
Fourth row: e"^e~7, invariant mass of the two electrons on different horizontal and vertical scales. 
Symbols as in Fig. [TJ 
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2540607-003 




p(J/\|/) (GeV/c) mass(GeV) cos0(£+) cos0 



FIG. 3: Top to bottom: 77, Stt^, tt^tt~tt^, tt^tt^j, e~^e~'y. For each channel, left to right: J/ip 
momentum, r] mass, polar angle of the positive lepton from the J/tp decay, and polar angle of an 
r] decay product (most energetic shower for 77 and Sir^, positive track for 7r"^7r~7r'^, 7r"'"7r~7, and 
e+e~7). Symbols as in Fig. [TJ 
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FIG. 4: Comparison of the results obtained in this analysis with the most precise measurements 
from other experiments 0, Q (top), and the PDG 2006 global fits 0]. 
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